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A New Method for Characterizing Very-Low-Mass Companions with Low 
Resolution Near-Infrared Spectroscopy 

Emily L. Rice^’^, Rebecca Oppenheimer^, Neil Zimmerman^, Lewis C. Roberts, Jr.^, Sasha 

Hinkley^’® 

ABSTRACT 

We present a new and computationally efficient method for characterizing very low 
mass companions using low resolution {R ~30) near-infrared {YJH) spectra from high 
contrast imaging campaigns with integral field spectrograph (IFS) units. We conduct 
a detailed quantitative comparison of the efficacy of this method through tests on sim¬ 
ulated data comparable in spectral coverage and resolution to the currently operating 
direct imaging systems around the world. In particular, we simulate Project 1640 data 
as an example of the use, accuracy, and precision of this technique. We present results 
from comparing simulated spectra of M, L, and T dwarfs with a large and hnely-sampled 
grid of synthetic spectra using Markov Chain Monte Carlo techniques. We determine 
the precision and accuracy of effective temperature and surface gravity inferred from 
fits to PHOENIX dusty and cond, which we find reproduce the low-resolution spectra of 
all objects within the adopted flux uncertainties. Uncertainties in effective temperature 
decrease from ±100-500 K for M dwarfs to as small as ±30 K for some L and T spec¬ 
tral types. Surface gravity is constrained to within 0.2-0.4 dex for mid-L through T 
dwarfs, but uncertainties are as large as 1.0 dex or more for M dwarfs. Results for ef¬ 
fective temperature from low-resolution Y JH spectra generally match predictions from 
published spectral type-temperature relationships except for L-T transition objects and 
young objects. Single-band spectra (i.e., narrower wavelength coverage) result in larger 
uncertainties and often discrepant results, suggesting that high contrast IFS observing 
campaigns can compensate for low spectral resolution by expanding the wavelength 
coverage for reliable characterization of detected companions. We hnd that SNR~10 is 
sufficient to characterize temperature and gravity as well as possible given the model 
grid. Most relevant for direct-imaging campaigns targeting young primary stars is our 
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finding that low-resolution near-infrared spectra of known young objects, compared to 
field objects of the same spectral type, result in similar best fit surface gravities but 
lower effective temperatures, highlighting the need for better observational and the¬ 
oretical understanding of the entangled effects of temperature, gravity, and dust on 
near-infrared spectra in cool low gravity atmospheres. 

Subject headings: infrared: stars — stars: late-type — brown dwarfs — techniques: 
spectroscopic 


1 . Introduction 

Extra-solar planets have been indirectly detected for over 20 years, and there are more nearly 
2,000 known systems, with over twice as many Kepler candidates awaiting confirmationli] Indirect 
detection techniques can provide mass, radius, and bulk composition measurements for extra-solar 


imag i ng of extra-solar p lanets (e.g., 


2010 

2007 


Currie et ah 


Knutson et ah 


Fortnev et al 

2008: 

Soiegel & Burrows 

2012) 

Kalas et al. 

2008: 

Vlarois et al. 

2008, 

2010: 


Lagrange et al 


2014|j a n d studies of transiting planets in secondary eclips e le.g.. iDeming et al 


20081 : ISing Sz Lopez-Moralesl l2009l : iFortnev et ahl 120131 1 have expanded the 


field to include pseudo-spectra l studies via photom etric observations made at multiple wavelengths, 
including narrow band filters (l.lanson et al.ll20iril ). 


Integral field spectroscopy will rapidly expand the study of substellar companion and exoplane- 


atmospheric properties, even within the same planetary system ( 

Patience et al. 

2010: 

Bowler et al. 

2010: 

Barman et al. 

2011a: 

Konopackv et al. 

2013: 

Oppenheimer et al. 

2013: Chilcote et al. 

20151. 


Several ground-based high-contrast imaging instruments have been specifically designed to directly 
image extra-solar planet s, including: the Near-Infrared Coronagraphic Imager (NICI) on Gemini- 
South ( Chun et al.ll2008l ). the High Contrast Instrument for the Subaru Next Generation Adaptive 
Optics (HiCIAO) and later the Coronagraphic High Angular R esolution Imaging Spectrograph on 
(CHARIS) on Subaru ( Suzuki et al. 201C; McElwain et al. 2nidL Project 1640 (P1640~) o n the Palo- 
mar Observatory Hale 200-inch telescope ( Hinklev et al. 2011b : Oopenheimer et al.ll2012 ). the Gem¬ 


ini Planet Imager (CPI) on Gemini-South (iMacintosh et aLll2014lL and the S pectro-Polarimetric 


High-contrast Exoplanet Research (SPHERE) on the VLT (jBeuzit et al. 2008) 


P1640, CPI, and SPHERE include integral field spectrograph (IFS) units that produce low 
resolution spectra (spectral sampling R=[A/AA]~30-60) of directly imaged companions in the near- 
infrared (~l-2 pm). The high contrast ratios achieved by these instruments greatly expands the 


^http://exoplanetar chive, ipac.caltech.edu/ 
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parameter space of the companions for which spectro scopy can be obtained, i ncluding M, L, &: T 


dwarf and planetary companions around BAFG stars ( 

Zimmerman et al. 

2010: 1 

[inklev et al. 2010. 

2011a: Roberts et al. 2012: Hinklev et al. 2013: 

Chilcote et al. 

2015: 

CreoD et al 

20151. and simul- 

taneous observations of multiple-planet systems 

(OoDenheimer et al. 

2011 

1), all at Solar System-like 


separations (4-40 AU). The low resolution of spectra obtained by these instruments are crucial for 
confirming the nature of the companions and for constraining their atmospheric properties and 
formation scenarios. 


The primary purpose of IFS units on high contrast instruments is to distinguish faint compan¬ 
ions from residual l ight from the primary star tha t has been diffracted through the instrument in 
a complex pattern (jOpDenheimer fc Hinklevi 120091 1 . The IFS exploits the wavelength dependence 
of the diffraction in order to separate speckles from faint objects, which will rema -in at the same 
posit ion on the image plane while the speckles change position with wavelength (jSparks fc Ford 
2002I 1. The spectral sampling and coverage of IFS units is motivated by speckle suppression rather 
than spectral characterization of detected companions and limited by the engineering complexities 
of the lenslet arrays and the pixel and chip size of infrared detectors. 

Very low-mass stars, brown dwarfs, and planetary-mass companions have potential to be effi¬ 
ciently characterized by comparing observed low resolution spectra to synthetic spectra from model 
atmospheres, as long as the selected wavelength regime contains significantly broad spectral fea¬ 
tures that are sensitive to atmospheric properties like effective temperature, surface gravity, and 
metallicity. Spectral fitting thus avoids assumptions about age, ra d ius, and mass, as well as th e 
evolutionary i nodels that relate these quan tities, e.g.. iBurrows et al.l (jl997l ): IChabrier et al.l ([20001); 
Baraffe et al.l ([20021 1: iFortnev et al.l ([20071 1. As higher-contrast and lower-mass objects are imaged 
and observed spectroscopically, the ability to conduct such comparisons efficiently and reliably is 
needed. 


We have developed a method to derive atmospheric parameters and uncertainties by fully 
mapping a broad temperature-gravity parameter space using thousands of synthetic spectra and 
linearly interpolating between calculated models. We compare simulated P1640 spectra (described 
in Section [2[) of field M, L, and T dwarfs and young M and L spectral type objects to synthetic 
spectra calculated with the PHOENIX model atmosphere code (Section [3|) using a Markov Chain 
Monte Carlo (MCMC)-based fitting procedure described in Section 01 Section [5] presents the re¬ 
sulting probability distribution functions and best fit model spectra, describes the constraints on 
parameters of effective temperature and surface gravity, compares results obtained using different 
subsets of the simulated spectra, and discusses the dependence of best fit parameters and uncer¬ 
tainties on the noise level in the simulated spectra. The conclusions of this work are presented in 
Section [H 
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2. Spectral Templates and Simulated Spectra 

In order to test our fitting procedure over a range of effective temperatures (hence companion 
masses), we constructed low resolution (AA~300 A) spectra from moderate-resolution, high signal- 
to-noise ratio (SNR) observed spectra of field M, L, and T dwarfs and young M and L dwarfs. The 
original observational spectra are henceforth referred to as “template” spectra and the constructed 
lower-resolution spectra as “simulated” spectra. The creation of simulated low-resolution near- 
infrared spectra from the template spectra is described below. 


2.1. Spectral Templates 


The sample of objects used as spectral templates is summarized in Table [T] Their spectral 
types range from Ml to T4.5, spanning the range of effective temperatures expected for low mass 
compan ions that can be detec ted by the current generation of high contrast integral field spectro¬ 
graphs (jBeichman et al.ll2010l f. All template spectra were obtain ed with the SpeX spectrometer on 
the NASA Infrared Telescope Facility fIRTF. lRavner et al.ll2003l f in cross-dispersed mode, covering 
0.81-2.4 /un with a resolving power {R=\/AX) of at least 2000, with the exce ption of the spectr um 
of PSO J318.5338—22.8603, which was obtained on the GNIRS spectrograph ( Elias et al. 2006) on 
the Gemini-North 8.1-m Telescope in cross-dispersed mode, covering 0.95-2.5 /rm with R ~1700. 

High-contrast direct-imaging surveys will primarily target young stars for which low-mass 
companions are likely to be still contracting, thus they have lower surface gravity than field-age 
objects of the same mass. Template objects include both field (ages >1 Gyr) M, L, and T dwarfs 
and you ng (^10-100 Myr) M and L dwarfs. Spe ctra of field objects are from the IRT F Spectral 


2010cl ). an LO7 ( Kirkpatrick et al.ll2006l ). an Lhy ( Fahertv e 


Librar y (IGushing et al.ll2005l: iR.avner et al.M2009l ). The young objects are a n M8.5o' (IRice et al 


complete description of the 7 gravity suffix), and an L7 VL-G (jAllers fc Liu 


al 


2013. see 


Gruz et al. 


2013 


200s for 


Liu et al 


20131). 


These young objects are high probability memb e rs of 0 Pictoris, Tu c ana-Horologium, AB Doradus 


and (3 Pictoris, respectively (jRice et al.l l2010d : iGagne et al.l l2014l : iFahertv et al.l 120131 : iLiu et al 


2OI3I I and are among the closest photometric and spectroscop ic analogs of direct l y-imaged gas 


giant planets for which detailed analysis is possible (see, e.g., 
2013l i. 


Fahertv et al.l 120131 : lAllers &: Liu 


2.2. Simulated High Contrast IFS Spectra 

We simulate low-resolution near-infrared spectra of very-low-mass stellar and substellar com¬ 
panions that can be obtained by high contrast IFS units in order to test the accuracy and precision 
of physical parameters inferred from model fits to these spectra. We report results for spectra rep¬ 
resentative of the resolution and wavelength range for Project 1640, which has a single observing 
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mode covering the Y JH band with a sampling of AA rsj 300 A. 


Project 1640 has operated since 2008 on the Hale telescope in two distinct phases, Phase 1 
and Phase 2, separated by a substa ntial instrument upgrade (new detector, filters, and increased 
sensitivity. lOppenheimer et alJl2012l ). Phase 1 (July 2008 to June 2012 ) produced JH (1.1-1.8 um ) 
spectra at 300 A pixel“^ resolution with one instrument setting (see iHinklev et al.ll2008l . l2011bl ). 
Phase 2 upgrades (begun June 2012) extended the wavelength coverage to YJH (0.9-1 . 8 jurn ) 
and slightly increased the spectral sampling to 286 A pixel” ^ ( Oppenheimer et al.l 120121 . l2013l ). 
Simulated P1640 spectra were created from SpeX/IRTF template spectra for both phases, but the 
results of this paper focus on fits to Phase 2 spectra. 

The first step in the creation of simulated P1640 spectra is to define the wavelength channels 
based on the data reduction and spectral extractio n procedures developed by the instrument team 
(IZimmerman et al.ll201ll : lOppenheimer et al.ll2013l ). At each wavelength channel, the flux density 
{F\) value of the simulated spectrum was set to the median of the flux density from the template 
spectrum within the wavelength range of the channel. The quoted resolution of the binned spectra 
represents the spectral sampling rather than the full-width at half-maximum of the instrumental 
profile. The simulated spectra were then normalized to the maximum value in the P1640 wavelength 
range. The template IRTF spectra and simulated P1640 spectra are plotted in Figures [T] and [2l 

The simulated spectra are significantly higher SNR than the typical spectra of very-low- 
mass companions obtained with ground-base d high contrast integral field spectrographs (e.g., 
Oppenheimer et al.ll2013l : IChilcote et al.ll2015l ). Therefore, for the primary analysis described in 
Section m the uncertainty on the simulated spectra was set to SNR=10 at the peak flux value and 
assigned a constant noise at each flux point, resulting in lower SNR for lower flux values. The ac¬ 
tual uncertainties on these flux points (i.e., errors on the template spectra flux added in quadrature 
within each wavelength channel) the are much smaller; therefore, we test the dependence of best 
fit parameters and uncertainties on SNR via several “resampling” tests described in Section 15.71 

The simu lated spectra do not explicitly in clude any instrumental systematics such as speckles or 
cross-talk f see lOppenheimer &: Hinklevll2009l ). The effects of residual uncorrected wavefront aberra¬ 
tions, or speckles, coinciding with a detected companion in the extracted image cube is prohibitively 
time-consuming to predict and incorporate into the our anal ysis, although they are minimized with 


current data reduction a nd spectral extraction procedures (jPuevo et al.l 120121 : lOppenheimer et al 
20131 : iFergus et al.ll2ni4l ). Another potentially significant instrumental effect is cross-talk, or light 
from different wavelengths leaking in to neighboring spectra on the detector. For Phase 2 data, 
cross-talk has been reduced to a maximum of 1.5% and average of 0.4% across all wavelengths with 
better focusing of the inst rument, an improved blocking filter, and optimized extraction methods 
(jOppenheimer et al.ll2013l ). 













































3. Model Atmospheres and Synthetic Spectra 


Model atmospheres and synthetic spectra were calculated using the PHOENIX atmosphere code 
(jHauschildt et al.lll997l . Il999l l. For an effective temperature range 1400-4500 K, appropriate for M 
and L dwarfs, the dus ty ver sion of the PHOENIX code was used, as presented in iRice et alJ (j2010al h 
Schlieder et al.l (|2012l ). and [Roberts et al.l (j2012l ). The cond version of the PHOENIX code was used 
for an effective temperature ran ge of 950-2000 K, appropriate for mid to late L dwarfs a nd early 
to mid T dwarfs as presented in iRice et alJ (l2010bl i and similar to iDel Burgo et alJ (120091 1 . 


The dusty and cond versions of the PHOENIX models represent limiting cases of dust treatments. 
In the dusty version, dust is created in thermo-chemical equilibrium, assumed to be made of spherical 
single-species grains with a power-law size distribution, and remains static in the atmosphere as a 
source of opacity. In the cond version, dust is created identically but the opacity is not considered 
in calculation of the emergent spectrum, thus mimicking perfectly efficient removal of the dust to 
layers below the photo sphere. These lim iting cases have been shown to work well f or objects hotter 


than mid-M type le.g.. iRice et al.ll2ninal l and cooler than the L-T transition (e.g.. lDel Burgo et al 


20091 : mice et al.l 1 201 Obi ) . While more sophisticated dust treatments have been implemented in 


other atmosphere models, we use PHOENIX models because we have access to a complete and finely 
sampled model grid over the required range of effective temperature and surface gravity. 


To test the veracity of our PHOENIX dusty and cond synthetic spectra compared to low-resolution 
near-infrared spectra of M, L, and T dwarfs, we added a nuisance parameter to the fitting procedure 
described in Section |4] below. This parameter is added in quadrature to the actual uncertainty on 
the binned flux points (i.e., not the adopted uncertainty for SNR<10) and allowed to vary. The 
Te// and log(g) results from these fits marginalized over the tolerance parameter were within the 
1-cT uncertainties of the reported results for all spectra and the nuisance parameter was always less 
than 0.10, with typical peaks around 5%. This indicates that the systematic differences between 
the models and the simulated spectra are ~5% and as little as 1-2% for the field M dwarfs. 


For the model grids, atmospheric temperature-pressure structures were calculated in intervals 
of 50 K and with surface gravities ranging from log(g)=3.0-6.0 [cgs] in 0.1 dex intervals, resulting 
in 1952 dusty models and 682 cond models. We include surface gravities up to log(g)=6.0 (some 
unphysical) in our model grid in order to test the inherent sensitivity of the spectra to surface 
gravity in the absence of independent age constraints and to avoid a possibly artificial upper limit 
at log(g)=5.5. This also enables more reliable uncertainty estimates at high surface gravities because 
the probability distributions are not restricted to the physically allowed values. All atmosphere 
structures were calculated at solar metallicity. 


The calculation of a converged atmospheric structure with the PHOENIX code also produces 
a synthetic spectrum at 4 A sampling from 10 A to 5 /im; therefore, the models oversample the 
simulated P1640 spectra by a factor of ~70. Synthetic spectra matching the wavelength range 
and spectral sampling of PI640 are created following the same method used for simulated spectra, 
described in Section 12.21 above. Binned synthetic spectra for the range of effective temperatures 
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and surface gravities used in the fitting procedure are shown in Figure [3l Within the model fitting 
procedure (described in Section H]) the calculated synthetic spectra are first linearly interpolated 
in flux to produce spectra with the desired temperature and gravity values, then binned to match 
the wavelength coverage and spectral sampling of the simulated PI640 spectra as described in 
Section 12.21 


4. Spectral Fitting Method 


We use a two-step spectral fitting procedure to derive best fit physical parameters and un¬ 
certainties. The fitting procedure uses a grid of model spectra that are ~10 times higher spec¬ 
tral sampling than the comparison (i.e., observed or simulated) spectrum and calculated at fine 
enough sampling in parameter space to allow reliable linear interpolation in flux (as described in 
Rice et al.ll2010al and below). The linear interpolation is done on the fly by the fitting code, ef¬ 
fectively mimicking a continuous grid of model spectra in the specified parameter space. In the 
current analysis we consider model parameters of effective temperature and surface gravity. Future 
work will include super-solar abundances exp ected for extra-solar planets created via core accretion 

(e 




Fortney et al.l 20081 : lOberg et a 


ofiRice et al.l pOlOal ) and updated in 

(|20ld v 


201111. The f itting procedure is derived from the methods 


Roberts et al.l ( 20121 1. iHinklev et al.l ( 201.4 ). and ICrepp et al 


The first step of the procedure uses the calculated model grids described in Section [3l For 
each synt hetic spectrum we ca lculate a goodness-of-fit value based on the weighted sum of squared 
errors (see lCushing et al.ll2008l l. This calculation provides a coarse overview of the entire parameter 
range. Example results of this step are presented in Figure |4] for the Y JH simulated spectra of the 
field M5, LI, and T4.5 objects. 


The second step of the procedure uses a Metropolis-Hastings-based MCMC method to derive 
the posterior probability distribution over the model parameters. We assume an uninformative 
prior, i.e., one that is uniform over the range of parameters in the model grid. The uniform prior 
ignores prior knowledge about the objects that would make certain parameter regions higher or 
lower probability. We use it in order to present the most broadly applicable results and uncertainties 
without incorporating potentially inaccurate prior assumptions. The chain is initialized using the 
parameters of the model with the lowest (best) goodness of fit values from the previous step. 
Subsequent models are generated on the fly via linear interpolation of flux points from the pre¬ 
calculated grid in order to analyze an effectively continuous grid of models. We tested the veracity 
of linear interpolation in temperature by comparing spectra created via linear interpolation between 
two models 100 K apart to the calculated model with the same temperature. The mean difference 
for each flux point was 0.7% for the dusty models and 1.4% for the cond models, each with (T~2% 
so we expect the systematic differences from linear interpolation to be at most 1-2%. 

The uncertainties on the simulated or observed spectra are used to calculate the for the 











































samples drawn from the posterior distribution. The “jump size” (width of the normal distribution 
from which the new model parameters are randomly selected) was adjusted for each object to 
produce an acceptance rate (number of accepted jumps divided by the total number of links) of 
~0.3-0.4. The jump sizes varied from 12.5 to 1450 K in effective temperature and 0.025 to 2.9 dex 
in surface gravity, with larger jump sizes for earlier spectral types. If a jump took the parameter 
outside the calculated model grid, then the interpolation step sets the parameter to just inside the 
grid. 

The MCMC chain length for each spectrum is 10® links, where each link is a sample from 
the posterior distribution, whether or not the new sample is accepted into the chain. Typically 
in MCMC analysis the first ~10% of the links are considered “burn-in” and not included in the 
final analysis; however, initializing the chain with the parameters with the lowest global goodness 
of fit values effectively removes the need for a burn in period. We confirm this by comparing the 
mean and standard deviations of the first 10® steps in the chain with the 9x10® post-burn links, 
which were nearly identical. Therefore the complete 10®-link chains are used in analysis presented 
in Section [5l 


For preliminary tests of the htting procedure, spectra were fit with the both dusty and the cond 
models with the appropriate temperature ranges for each dust treatment (described in Section [3]). 
This, as expected, resulted in poor fits to the M dwarfs using the cond models and to the T dwarfs 
using the dusty models. For the L dwarfs, the dusty model fits resulted in qualitatively better fits 
to the data; therefore, only the dusty models were used for the complete analysis of the L dwarfs. 
Preliminary tests were also run with fewer links (10^-10®) with similar best ht results, indicating 
that the 10®-link chains are converged. 


P1640 spectra include the H 2 O absorption band at ~1.4 /rm that is present in objects ~M4 and 
later as well as in the Earth’s atmosphere. Because the signal is low and the correction for telluric 
absorption is often poor, this region is typica lly excluded from analysis for observed spectra (e.g.. 


Roberts et al. 


2012 ; 


Oppenheimer et al.l 120131 ). Therefore, we also fit “trimmed” simulated P1640 
spectra by excluding four flux points in the ~1.4 fim H 2 O absorption band. GPI and SPHERE 
cover the YJHK spec tral range with several individual hlters (e.g., Y, J, H, Kl, and K2 for GPI, 


Macintosh et al.ll2008h . Therefore we ht the simulated P1640 spectra as separate TJ-band and 


Fl-band spectra. Results of these hts are discussed in Section [5l3] below. Future papers will present 
this analysis for the exact wavelength coverage and spectral sampling of GPI and SPHERE. 


5. Results 

5.1. Probability Distributions 

Results from the spectral htting routine described above are probability distributions with 
10® values for each parameter. The probability distribution for hts to the simulated spectra from 
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the field M5 spectral template are shown in Figure O for the field LI in Figure [H for the field 
L8 in Figure [71 and for the field T4.5 in Figure [H In all figures the Y JH results are shown in 
purple, the YJ results in blue, and the H results in red. The lower-left panel shows 1-, 2-, and 
3-cj contours of the posterior distribution functions for both parameters, while the top and right 
panels shows the histograms for temperature and surface gravity, respectively, marginalized over 
the other parameter. 

The probability distributions for the field M5 are presented in Figure [5l The Y JH and 

Y J-band spectra produce similar temperature results, but the 77-band distribution has broad, 
flat peak with a long tail toward higher temperatures. None of the distributions shows a clear 
peak in surface gravity. None of the distributions shows a clear peak in surface gravity, with the 

Y JH distribution increasing toward unphysically high surface gravities and the Y J- and i7-band 
distributions increasing slightly at values of log(g)<4.0 dex (cgs). 

The probability distributions for the field LI (Figure [6]) are qualitatively similar to the M5 in 
temperature (~600 K cooler), with the YJH and TJ-band spectra producing nearly overlapping 
symmetric and clearly peaked distributions, but the J7-band distribution has a broad peak at hotter 
temperatures and low plateau at lower temperatures compared to the YJH and YJ peaks. In 
surface gravity, the Y JH and J7-band distributions show clear peaks, while the J-band distribution 
is fiat. Both the YJH and J7-band spectra show sharp decreases in probability at surface gravities 
>5.7, which are unphysically high values but allowed in the model grid. 

The probability distributions for the field L8 are presented in Figure [71 The distribution for 
the Y JH spectrum is narrow and symmetric in both temperature and gravity. The Y Jdistribution 
has a broader, but still single-peaked, distribution in temperature but increases toward unphysical 
high values in surface temperature. The J7-band distribution is broad in temperature, covering 
the entire range of both the YJH and YJ distributions in temperature. In gravity the J7-band 
distribution is consistent with the YJH distribution with a long, flat tail toward lower values. 

The distributions for the T4.5 object (Figure [8]) are the smoother and more symmetric in 
both parameters than those for the earlier spectral types. The peaks in the effective temperature 
distributions for each spectrum are offset by ~50 K but are within the 1-a uncertainties of the 
other distributions. The surface gravity distributions are similarly offset, J7-band to higher values 
and Y J-band to lower values, but are also consistent within the 1-a uncertainties. Both parameters 
are better constrained by complete YJH spectra, with temperature more constrained by the YJ 
fit and gravity more constrained by the J7-band ht. 

Comparing Figures [SHU h becomes apparent that the low-resolution near-infrared spectra of 
very low mass objects are increasingly sensitive to temperature at later spectral types. Furthermore, 
the distributions from fits to the Y J-band are generally more consistent with those of the complete 

Y JH spectrum. The YJ spectrum distribution for the L8 template object is an exception, but 
this is also the temperature regime where the dusty PHOENIX models are not expected to accurately 
reproduce cool atmospheres. In surface gravity, the probability distributions also show tighter 
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constraints for later spectral types, although these distributions sometimes extend into unphysical 
regimes of log(g)>5.6. 


5.2. Best Fit Spectra 

We adopt the best fit parameters of the 50% quantile value for each parameter distribution, 
which corresponds to the peak of the histograms only for symmetric distributions. The uncertainties 
on the best fit parameters are defined by the 16% and 84% quantiles of the distributions, which 
correspond to 1-cr width for a purely Gaussian posterior. Results are listed in Table[2]and discussed 
below. 

The effective temperature and surface gravity values of the best fit model spectra for each 
simulated P1640 spectrum are listed in Table [21 and the corresponding synthetic spectra are over¬ 
plotted with simulated spectra in Figures [9] and [101 Best fit parameters were determined via the 
two-step method described in Section [H above for four versions of the simulated spectra: complete 
YJH spectra, “trimmed” YJH spectra, and separate YJ and H spectra. The best fit results for 
the complete and trimmed YJH spectra are effectively the same (see Section 15.31 below) so only 
the complete YJH results are reported and plotted, as solid purple lines, in Figures P and [TOl Blue 
dashed lines represent the best fit synthetic spectrum with model parameters determined using 
only the FJ-band flux points (0.995-1.31 ^m) and red dot-dashed lines, the H-hand flux points 
(1.45-1.80 /rm). All three best fit spectra are plotted over the complete YJH wavelength range of 
the simulated spectra. 

The best fit synthetic spectra from the YJH fits reproduce the overall spectral shape and 
broad absorption features, within the assigned error bars of the simulated spectra, for all spectral 
types. The only marginal exceptions are 1-2 points in the the trough of the 1.4 /un H 2 O absorption 
band for the field L5 and L 8 objects, one point in the F-band peak of the T4.5 dwarf, as well as 
two points around 1.4/im and the slope on either side of the //-band peak for the LO 7 object. This 
is remarkable especially considering the limiting-case dust treatments of the PHOENIX dusty and 
cond models and the discrepancy of the model parameters with previous measurements from the 
literature (Sections 15.41 and 15.5|) . Therefore we advise that even good fits should be treated with 
caution when interpreting fits to low-resolution near-infrared spectra for very low mass objects. 

The fits to different subsets of the flux points, based on traditional filters separated by the 
telluric absorption band at 1.4 /un, test the consistency of the best ht parameters from these 
different bands. For the M dwarfs, all of the hts produce similar results, with the //-band fit 
slightly over predicting the hrst two or three F-band flux points for the M5 and M 8 . For the 
L dwarfs, the fits become less consistent, with the //-band parameters over predicting the flux for 
the F J-band points and the F J-band ht under predicting the hux at //-band, especially for the L 8 
object. The T dwarfs are slightly better, with the only problems being the F J hts under-predicting 
the red end of the //-band for the T2 and slightly over predicting the hux for the T4.5, and the 
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-ff-band fit slightly under-predicting the Y- and J-band peaks for the T4.5. 

The offsets between the synthetic spectra from Y J- or H-hand hts and the rest of the simulated 
spectra are more signihcant for the young objects, particularly for the low gravity L spectral types. 
Figure [IHl shows the best ht spectra for the young objects (left column) with the best ht spectra 
for held objects of similar spectral types (right column). For the hotter objects (MS.hy and LO 7 ) 
the yj-band hts only slightly under-predict the Ff-band hux points (blue dashed line), while the 
FF-band hts are more signihcantly different, especially for the LO 7 spectrum (red dot-dashed line). 
The single-band best ht spectra for the cooler objects (L 57 and L7 vl-g) always under-predict the 
hux in the other band, most substantially for the TJ-band ht in the H band (blue dashed line). 

Overall, the low-resolution hts to the complete YJH spectra reproduce the simulated low 
resolution near-infrared spectra for all spectral types, despite using atmosphere models with limiting 
cases of dust treatment to produce the synthetic spectra. Fits to an individual bands result in 
signihcant discrepancies with hux points in the other band, suggesting that a broad wavelength 
coverage is best for reliable characterization of low mass companions via low resolution spectra. 


5.3. Effect of Trimming Spectra 


The H 2 O absorption band at ~1.4 /rm is typically low SNR and most likely to be affected by 
systematic errors induced by telluric calibration, which are difficult to quanti fy. These wavelength 


points are typically exclu ded from analysis for observed P1640 spectra (e.g.. lEoberts et al.l 120121 : 


Oppenheimer et al.ll2ni3l ). Therefore, we ht “trimmed” simulated P1640 spectra by excluding four 


hux points in the ~1.4 /xm H 2 O absorption band to determine the effect of removing these hux 
values on the best ht model parameters. 


The resulting best ht parameters were identical within the uncertainties in both cases, with 
the largest differences being 15 K for the M5 object and 0.15-0.25 dex for the held L5 and L 8 
objects and the young M 8.57 object. The temperature difference is 10 times smaller than the 1-a 
uncertainty dehned by the 16% and 84% quantile values. The gravity differences are comparable 
to the 1-a uncertainties for those objects. Uncertainties derived from the probability distributions 
were also very similar, within 15 K and 0.04 dex log(g), except for the log(g) of the M 8.57 object, 
with the 1-a quantile extending 1.4 dex lower in surface gravity for the trimmed spectrum. 

These results suggest that excluding low SNR hux measurements in the water band have little 
effect on the parameters of very low mass objects inferred from model hts to low resolution near- 
infrared spectra, but they might be marginally signihcant for hotter low-gravity objects. Results 
for the young objects are discussed in more detail in Section [5ffil below. 
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5.4. Effective Temperature 


Figure [TT] shows the effective temperatures of the best fit synthetic spectra as a function of 
spectral type for the field objects (filled symbols) and young objects (open symbols) for each of the 
three versions of the spectra: V J (upward blue triangle), H (downward red triangle), and complete 
YJH (purple square). The triangle points are offset slightly in spectral t ype for clarity. The gray 
shaded region sho ws temperatures deriv ed for the same spectral types from Luhman fc Rieke (jlQhd ) 
for M dwarfs and [Stephens et al.l (j2009l l for L and T dwarfs with a range of ±100 K. 


For model fits to low-resolution near-infrared spectra to be considered reliable, they should 
be consistent with the physical parameters derived from other methods. Figure [TT] shows that the 
best fit results from low-resolution near-infrared spectra (symbols) match the empirically derived 
temperatures (gray shaded regions) to within the uncertainties for the field M dwarfs and the T4.5 
dwarf, but the model fits result in hotter temperature for the field L5, L 8 , and T2 dwarfs, by as 
much as 500 K (for the L 8 dwarf). This is the temperature regime in which dust is expected to 
condense from the photosphere, a dynamic and complex process that is not modeled in detail with 
the limiting cases of dust treatment in the PHOENIX dusty and cond model atmospheres. 


The best fit effective temperatures from different segments of the spectra are generally consis¬ 
tent for a given object within the 1-a uncertainties, with notable exceptions of the Y J fits for the 
L7 VL-G, L 8 , and T2 objects. The TT-band fits produce the highest temperatures for almost all of 
the field objects and lower temperatures for the young L dwarfs. 


Uncertainties in effective temperature are large (up to ±500 K) for the earliest spectral types 
considered in our analysis, but these are still 10-20% precision. Statistical uncertainties are as 
small as ±30-50 K for the mid to late L dwarfs (field and young), although these values are also 
systematically hotter than the empirically derived temperatures, particularly for the field objects. 
As might be expected, the uncertainties are smallest for the complete YJH spectra, but in some 
cases the uncertainties on individual bands are nearly as small, particularly for the Y J-band of the 
M 8 . 57 , LO 7 , L5, and L7 objects and for the F7-band of the L 57 , L7 VL-G, and T2 objects. 


The best fit effective temperature for known young objects are shown on Figure [TT] as open 
symbols. All young object spectral fits produce cooler temperature best fits than corresponding fits 
to spectra of field objects with similar spectral types. This is consistent with the idea that enhanced 
dust in cool, low-gravity atmospheres causes the near-infrared spectrum t o appear more red (i.e., 
cooler) than would be predicted by the optically - defined spectral typ e fe.g.. iMetchev fc Hillenbrand 


20061 : iLuhman et al.ll2nr)7l : [Barman et al.ll2011bl : iBowler et al.ll2013l ). Indeed, direct empirical com¬ 


I 


parison of bolometric luminosities for young and field objects of the same spectral types suggest 
that, because they are similar but the young objects should have inflated radii, the young objects 
must have lower effective temperatures than their field-age spectral type counterparts (Filippazzo 
et ah, in prep.). 
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5.5. Surface Gravity 


Substellar-mass objects never reach a stable main sequence; instead they cool, fade, and shrink 
for their entire lifetimes. Effective temperature and luminosity are degenerate with age and mass, 
but surface gravity is uniquely and significantly low (<5.0 dex cgs) for young objects (<100s Myr). 
We include four known young substellar objects in our sample to test whether model fits are 
sensitive enough to surface gravity to distinguish young from old objects, or even to constrain age 
via comparison with evolutionary models, using low resolution near-infrared spectra. 


Figure [12] shows the surface gravity of the best fit model as a function of spectral type with the 
same symbols as in Figure [11] described in Section [5.41 The shaded regions represent predictions for 
field (age> 5 00 M yr, solid) and young (age 5-120 Myr, ha tched) objects from e v olutio nary models: 

( 2000 1 for Tpff> 3000 K, Baraffe et al. (2002) and Chabrier et al. (2000) [DUSTYOO] 


Siess et al. 


for 3000>Te//>1500 K, and iBaraffe et al.l ([2003]) [CONDO] for Te//<1500 K. The hatched gray 
regions shows predictions from the DUSTYOO models for effective temperatures ~2500-1500 K, 
which approximately corresponds to the spectral types of young objects in our sample. 


Broadly, evolutionary models predict the surface gravity for field early-M dwarfs to be in the 
lowest part of the shaded parameter range on Figure [TJ] i.e., log(g)=4.9 dex. The oldest late- 
L dwarfs are predicted to peak at log(g)=5.4 dex, and the oldest, most massive T dwarfs reach 
log(g)=5.5 dex. For substellar-mass objects, younger objects will have lower surface gravities as a 
result of their larger radii, but objects older than 500 Myr will have log(g)>4.9 for the entire range 
of temperatures considered, according to both the DUSTYOO and COND03 evolutionary models. 


The best fit surface gravities are consistent with the values predicted by evolutionary models 
for field objects of all spectral types, within the relatively large 1-a uncertainties (from ±0.2 to 
±0.9 dex). For the M and early-L dwarfs, the YJH best fit surface gravity is closer to the predicted 
value than the surface gravities from individual Y J and ff-band fits. For field L and T dwarfs, the 
best fit surface gravity is within the range of values predicted by the evolutionary models, except 
for the individual Y J- and ff-bands for the LI object. This is likely a result of the broad spectral 
features being increasingly shaped by surface gravity, in addition to temperature, for the coolest 
objects. It is perhaps surprising that such low resolution spectra are sensitive to surface gravity 
at all, but it has been shown that the relative strengths of broad band absorption like H 2 O and 


collisionally-induced absorption from H 2 depend on surface gravity (jRice et al 


201 ll . Figure 6). 


The results from fits using different wavelength ranges for field objects are generally self- 
consistent to within the 1-cr uncertainties. The largest discrepancy is for the field L5, where the 
individual YJ- and iL-band fits indicate lower surface gravities (~5.0 dex cgs) and the complete 
YJH fit indicates a higher surface gravity (5.35 dex cgs). The three fits indicate similar effective 
temperatures (~2100“2230 K), with the individual bands at the hotter end of the range. The 
results are all within the l-cr uncertainties for each parameter so are not likely to indicate any 
systematic trend. 
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The results for best fit surface gravity for the young objects are more discrepant from pre¬ 
dictions of the evolutionary models and less consistent for fits from individual wavelength ranges. 
They are discussed in more detail in Section 15.61 below. 


5.6. Young Objects 


The primary targets of exoplanet direct-imaging campaigns are nearby young stars. Planetary- 
mass companions to young stars are more luminous than older objects because they are still 
contracting and radiating gravitational potential energy. Young brown dwarfs are similarly low 
mass and low surface gravity, but as free-floating objects they are amenable to a broader range 
of observations, both in wavelength coverage and spectral resolution. Therefore, young brown 


dwarfs are potentially i mportant analogs to directly-imaged exoplanets (see, e.g., iFahertv et al 


201J; lAUersJUiul 12013!). 


We include four conhrmed young brown dwarfs in our analysis in order to test whether model 
spectra reproduce the low resolution near infrared spectra of cool, low gravity objects and to 
compare the accuracy and consistency of best fit results for young and held objects with similar 
spectral types. The young template objects we use have opt ical or near-infrared sp ectral types of 
M 8 . 57 , LO 7 , L 57 /L 3 VL-G, and L7 VL-G (see Cruz et al. 2009 and Allers fc Liu 2013 for descriptions 
of the surface gravity suffixes) and medium resolution, high SNR spectra available in the literature 
(see Tabled] for references). 


Figure [l3] shows the posterior probability distribution functions for young M 8.57 and L 57 
objects and the held objects with the closest spectral types, with the YJH results in purple, the 
YJ results in blue, and the H results in red. The scales on each axis are identical in effective 
temperature and log surface gravity in order to show systematic difference between results for 
objects on the same optically-dehned spectral type but difference ages (held objects in top row and 
young objects in bottom row). For both young objects, the posterior distribution functions are 
centered on signihcantly cooler temperatures (ATe//=300-600 K) than for the comparable held 
object. For the young M 8.57 object, the YJH and YJ hts are consistent in temperature, but the 
i7-band ht is bimodal with the a signihcantly cooler component that is a very poor ht to the YJ 
spectrum. For the young LSy object, the complete YJH spectrum temperature distribution is 
consistent with the Y J and H distributions to within 1-a. 


Figure dH shows the posterior probability distributions for model hts to the L7 VL-G spectra. 
The distributions from hts to different spectral regions are substantially less consistent than for the 
held objects (Figures [SHH])- For example, the 2 -cj contours of the YJ-band (blue) and F7-band (red) 
distributions just barely overlap. The Ff-band distribution peaks at a surface gravity predicted for 
a young, late-type object (see Figure d2]), but the YJH and YJ-band distributions peak at higher 
surface gravities and have only a low tail toward lower values. This is further evidence that the 
PHOENIX dusty models do not accurately reproduce the observed spectra of young, late-type objects 
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Schmidt et al. 

2008: 

Mohantv et al. 

2010 : 

Patience et al. 

2012 ) 


Figure [TOl shows the simulated and best fit spectra for the young objects (left column) along 
with simulated and best fit spectra for field objects of similar spectral types (right column). The 
model spectra with the YJH best fit parameters are plotted as solid purple lines, the Y J-band 
fit as dashed blue lines, and the Ff-band fit as dot-dashed red lines. For all four young objects, 
the complete (and trimmed) best fit Y JH spectrum reproduces the simulated spectrum to within 
the adopted > 10 % uncertainties; however, the best fit parameters do not match predictions from 
evolutionary models, as discussed below. 


The field M 8 and young MS.Sy have similar optically-determined spectral types but the Y JH 
spectrum of the young object has a more steeply sloped Y J band, deeper H 2 O absorption, and a 
higher, more triangular FF-band peak. Fits to the complete YJH spectra are of similar quality, 
but fits to different wavelength ranges are less consistent for the young object than for its field-age 
counterpart. For the young object, the Y J-band fit slightly under-predicts the flux throughout 
the FF-band, and even the FF-band fit is slightly lower than the object’s flux beyond 1.7 //m. The 
FF-band fit over predicts the T J-band flux even more, well outside of the adopted flux uncertainties. 
This is qualitatively similar to single-band results for the field L5 and L 8 dwarfs, which are likely 
the dustiest of the field objects. 


The effective temperatures of the M 8.57 best fit spectra are 400-600 K cooler than those for 
the field M 8 (Figure [TT|). The best fit surface gravity for the young M 8.57 is comparable to that of 
the and field object in YJ- and FF-bands and match predictions from the DUSTYOO evolutionary 
models for 5-120 Myr objects (gray hatched region on Figure [12]). The best fit surface gravity 
for the YJH best fit spectrum is 1.0 dex (cgs) more than that of the field object, and is in fact 
higher than the maximum value predicted by evolutionary models for any substellar object (e.g.. 


Chabrier et al. 


2 OOOI : 


Baraffe et al 


2002 , 


2003). 


The complete Y JH fits to the simulated spectra of young L spectral types generally reproduce 
the flux points within the 1-cr error bars, except for the longest-wavelength FF-band point for three 
objects and several J- and FF-band points for the LO 7 object. For the LO 7 object, the YJH fit 
appears overall too flat at the blue end of the FF band, which corresponds to a high best fit surface 
gravity. The inconsistency between fits for individual bands is even more substantial for the young 
L objects than for field L dwarfs. For all three objects, the individual YJ and FF-band fits under 
predict the flux in the other band. The worst FF-band fit in the Y J-band is for the LO 7 , while the 
worst Y J-band fit in the FF-band is for the L7 VL-G object, which is the reddest of the young L 
dwarfs. 


All fits for the young Ls produce lower effective temperatures than for the field objects of 
similar spectral type, especially for the YJH and FF-band fits (Figure fTTI) . The best fit surface 
gravities for the LO 7 are similar to that of the M 8.57 in that the best fit parameters from the Y JH 
fits are too high, but the best fit parameters from the YJ- and FF-band fits are comparable to 
predictions from the DUSTYOO evolutionary models for young objects. The surface gravity results 
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for the later-type young L dwarfs are all higher than predicted by the evolutionary models, but not 
unphysical for field objects, and the //-band fits are closest to the predictions from evolutionary 
models (Figure [l2]), where the spectral shape is expected to be particularly sensitive to surface 
grav ity because of de creased collisionally-induced absorption from H 2 relative to H 2 O opacity (see 


e.g. 


Rice et alJl201ll .. Figure 6). 


Inconsistencies between model fits from individual bands, evident in both the probability dis¬ 
tributions and in the best ht spectra, likely stem from the dust treatment in the dusty models, 
which do not incorporate enhance d dust content, par t icularly in small gra ins, expected for low 
surface gravity atmospheres (e.g., ISchmidt et al.l l2008l : iMarocco et al.N2014l l. These results sup¬ 
port the growing consensus that ad d itional opacity source s and more sophisticated dust and cloud 


treatments (e.g., IWitte et al. 


2009; 


Barman et al.l 1201 lal : iMorlev et al.l 120121 . l2014l : lAllardI l2014l ) 


are necessary for modeling the atmospheres of young brown dwarfs as well as of directly-imaged 
exoplanets. 

These young, low-mass objects in particular illustrate the necessity of comparing observations 
in multiple bands in order to test the reliability of spectral characterization. Thus we caution 
against using results of model hts to low-resolution spectra to confirm or rule out youth of very 
low mass objects. S imilar results for the P1640 spectrum of k Andromedae B are described in 
Hinklev et al. ( 201,1 ). 


5.7. Signal to Noise Ratio 


The probability distributions, best fit parameters, and their uncertainties described above were 
calculated using the >10% uncertainties assigned to the simulated spectra described in Section [2]2j 
The noise level was selected such that the highest flux point of each simulated spectrum had 
SNR=10, which is e xpected to be typical for s pectra of exoplanets from high contrast integral held 


spectrographs fe.g.. lOpDenheimer et al.ll2013l : IChilcote et al.ll2015fl . The actual uncertainty on the 


binned hux points (i.e., the uncertainty on the observed hux points within the spectral range of 
each P1640 wavelength channel, added in quadrature) are signihcantly lower than 10% because the 
original SpeX and GNIRS template spectra have signihcantly higher SNR and we have not adjusted 
the binned hux values at all. This begs the question of whether the results described above are 
truly representative of model hts to SNR<10 spectra. 


Therefore we use two Monte Carlo simulations to test the dependence of the best ht parameters 
and their uncertainties on the adopted SNR of the simulated spectrum. The hrst Monte Carlo test 
was applied to to all spectral types with the adopted (SNR<10) uncertainties, and the second test 
used the M8, L5, and T4.5 simulated spectra to compare results from a range of different adopted 
uncertainty values. The complete Y JH simulated spectra were used for both SNR tests, and both 
tests compare simulated spectra to the calculated model spectral grid (i.e., 50 K intervals in effective 
temperature and 0.1 dex intervals in surface gravity. Figure [3]), without the interpolation used in 
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the MCMC fitting. 


5.7.1. Resampling versus real SNR 

The first SNR test, performed for all spectral types, resamples the simulated spectra by adding 
a random number to the binned flux values from the template spectra. The random number is drawn 
from a Gaussian distribution the width of the adopted noise level (SNR<10), hereafter called a 
“resampled” simulated spectrum. For 10“^ resampled simulated spectra, we calculate the goodness- 
of-fit to each calculated model spectrum, define the minimum goodness-of-fit value as the best fit 
parameters, and calculate the mean and standard deviation of the 10 ^ best fit values. 

The mean best fit parameters for 10^ resampled spectra of each template object are consistent 
with the YJH results listed in Table [2] to within one standard deviation. The largest discrepancies 
are for the M5 field object in surface gravity (different by ~0.5 dex) and the LO 7 object in effec¬ 
tive temperature (different by ~60 K), both still within 1-cr uncertainties of the MCMC results. 
Interestingly, the LO 7 object also has the worst “best fit” YJH spectrum (see Figure [inD. The 
uncertainties were generaliy about the same or smailer (by up to a factor of 2 ) for the resampied 
spectrum fits, but this is possibiy a consequence of using oniy the model grid instead of allowing 
for linear interpolation between calculated models. 


5.7.2. Dependence of Uncertainty on SNR 

For the second SNR test, we use the simulated spectra for the M5, LI, and T4.5 field objects 
to test the dependence of best fit parameters on SNR. As with the first test, the fiux points are ali 
resampied from within a Gaussian distribution the width of the noise on that flux point, this time 
for 20 values of maximum SNR, ranging from SNR=2 to 100. Each template spectrum is resampled 
lO'^ times for each SNR value, and the best fit model has the minimum goodness-of-fit value when 
compared to the resampled spectrum. The mean and standard deviation of the best fit parameters 
from 10^ trials are plotted as a function of SNR in Figure [151 

The best fit parameters for the resampled simulated spectra approach Te//=2950 K and 
log(g)=5.4 for the M5, Te//=2300 K and log(g)=5.5 for the LI, and Te//=1200 K and log(g)=5.3 
for the T4.5. These values are consistent with the MCMC YJH results within the 1-cr uncer¬ 
tainties, but some are offset from the closet model grid point to the MCMC result, likely because 
of the fundamental difference between MCMC and purely Monte Carlo techniques and the use of 
interpolation versus the calculated model grid. 

By SNR>5, the mean temperatures for all three objects converge to within a few percent 
of the SNR=100 value, and the surface gravities converges to the SNR=100 values for SNR>10. 
The uncertainties in surface gravity remain large (~0.15 dex) for the M5 dwarf, likely because the 
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low-resolution near-infrared spectra of M dwarfs are not as sensitive to gravity as the spectra of 
cooler objects are (i.e., Figure [5|). 


It should be noted that these tests assume a Gaussian noise model for the simulated IFS 
spectra, which may not be the most appropriate noise model for low SNR data. However, the 
Gaussian model is reasonable for a long chain of error sources with different characteristics, as is 
the case for high contrast IFS data. We expect the effects of the assumed error distribution to 
be smaller than the effect of the instrumental and the model systematics; therefore, we leave the 
evaluation of different noise models for a future paper. W e also ignore instrumental systematic 
uncertainties mentioned in Sections 12.21 and 15.31 However, lOppenheimer et al.l (j2013l i tested the 
fidelity of their S4 spectral extraction method by injecting fake sources with a T4.5 spectrum into the 
data cubes obtained for HR 8799bcde and comparing the extracted spectra to the input spectrum. 
The average deviation from the input spectrum over all wavelengths was ~2-9%, except for the 
innermost planet HR 8799e, which was 15%. Therefore we expect the instrumental systematics to 
be below the SNR~10 level we adopt for the simulated spectra. It should also be noted that the 
errors bars on Figure [15] are not representative of the actual uncertainty in the best fit parameters 
at high SNR, just the distribution of best fit parameters for resampled spectra, which will always 
approach a delta function at high SNR even for poor quality fits. 


6 . Conclusions 

We have developed a robust spectral fitting method and tested it for low-resolution {R ~30 
to 60) near-infrared (YJH) simulated spectra created from higher resolution observed template 
spectra of MLT dwarfs. We used the fitting method to explore the sensitivity of these spectra to 
effective temperature and surface gravity using the limiting-case dust treatments of the dusty and 
cond PHOENIX atmosphere models for M/L dwarfs and T dwarfs, respectively. 

The PHOENIX dusty and cond models reproduce the simulated spectra of both field and young 
objects when the entire wavelength regime is used in the fitting procedure, despite the limiting 
case dust treatments in these models. Best fit spectra determined using only one band {YJ or 
H) are typically within the adopted uncertainties of the complete simulated spectra for the field 
dwarfs, with some exceptions for the L5, L8, and T dwarfs. However, the single-band best fits 
generally do not reproduce the spectrum in the other band for spectra of young objects, which 
leads us to recommend caution when using single bands of low-resolution spectra to characterize 
young, low-mass objects. 

Our results indicate that low resolution near-infrared spectra are sensitive to temperature with 
precisions as good as ±50-100 K for L and T dwarfs, while the constraints are much looser (±100- 
500) for M dwarfs. The low-resolution near-infrared spectra are increasingly sensitive to surface 
gravity at later spectral types, although with lower precision (up to ~1 dex for hotter objects) and 
significant offsets from predictions of evolutionary models for the young objects. 





- 19 - 


The complete and trimmed Y JH spectra produced similar results and uncertainties, indicating 
that losing a few flux points to poor sky subtraction, speckle suppression, or other reduction and 
extraction effects will not significantly affect the characterization of the companion. 

The best fit parameters and uncertainties can be considerably different for the YJ and H 
subsets of the complete YJH spectrum, especially for temperature. Surface gravity results are 
generally consistent within the (large) uncertainties, but are discrepant from predictions of evo¬ 
lutionary models (see below). These results suggest that high contrast IFS observing campaigns 
should characterize detected companions with as broad a wavelength coverage as possible. In future 
analysis we will expand our simulated spectra into the K band and test the sensitivity and accuracy 
of different filter combinations to model parameters. 


The temperatures from the Y J spectral fits are typically more consistent with the temperatures 
from the full Y JH fits and closer to temperatures derived for objects of the same spectral types 
using other methods with the notable exceptions of the young LSy and L7 VL-G objects, for which 
the i7-band hts are closer to complete YJH results (Fig. [TT]) . 


The best fit temperatures are most discrepant from literature values for the field L5 and 
L8 objects, likely the result of our use of the PHOENIX dusty models. The best fit gravities are 
generally consistent with evolutionary models for the held objects but are signihcantly discrepant 
for the young objects. The complete YJH spectra hts to young objects always produced high 
surface gravities. Ev en Ff-band spectr a , which have been shown to be gravity sensitive at higher 


spect ral resolutions (jLucas et al. 


2001 


Luhman et al.l l2005l : lAllers et al.l 120071 : iRice et al.l l2010d . 


boilh . only produced marginally consistent results for the young L dwarfs (see Figure [T^ . Thus 
we caution against using model hts to low-resolution near-infrared spectra to conhrm or rule out 
youth of very low mass companions. 


We advise that even good hts should be treated with caution when interpreting hts to low- 
resolution near-infrared spectra for very low mass objects. The results for simulated spectra from 
L-T transition objects and for young objects in particular highlight the need for testing the more 
sophisticated dust treatments now available in cool atmosphere models (e.g., iMorlev et al.ll2012l j 
and for comparing best ht parameters over a broad wavelength range. 


We determine that a minimum of SNR~5 is required to reliably constrain the temperature of 
a low-mass companion and that SNR>10 is ideal for constraining surface gravity. We caution that 
these results are based on tests using objects with spectra that are reasonably well reproduced by 
the PHOENIX dusty and cond models and may not be representative of results for all directly-imaged 
high contrast companions. 


Future work will expand this analysis to simulated spectra for GPI and SPHERE, which cover 
YJHK with several individual hlters, and to include varying metallicity and dust treatments in 
the atmosphere models we consider. 
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Table 1. IRTF/SpeX Spectra 


Spectral 

Object 

R.A. 

Decl. 

Ref. 

Type 

Name 

(h, m, s) 

(“//') 


Field Templates 

Ml 

HD 42581 

06 10 34.6 

-21 51 53 

1 

M5 

G1 866 ABC’" 

22 38 33.7 

-15 17 57 

1 

M 8 

G1 752B 

19 16 55.3 

+05 10 11 

1 , 2 

LI 

2MASS J1439+1929 

14 39 28.4 

+19 29 15 

1 

L5 

2MASS J1507-1627 

15 07 47.7 

-16 27 39 

1 

L 8 

DENIS J0255-4700 

02 55 03.6 

-47 00 51 

1 

T2 

SDSS J1254-0122 

12 54 53.9 

-01 22 47 

1 

T4.5 

2MASS J0559-1404 

05 59 19.1 

-14 04 49 

1 

Young Templates 

M8.57 

2MASS J0608-2753 

06 08 52.8 

-27 53 58 

3 

LO 7 

2MASS J0141-4633 

01 41 58.23 

-46 33 57.4 

4 

L 57 /L 3 VL-G 

2MASS J0355+1133 

03 55 23.4 

+11 33 44 

5, 6 

L7±1(IR) VL-G 

PSO J318.5338-22.8603 

21 14 08.026 

-22 51 35.84 

7 


Refe rences . — (1) R.avner et al. ( 2009h . (2) Cushing et al 



1 4') iKirkpatrick et al 


il'] iLiu et al.l (1201311 




(5) iFahertv et al 


120051), (pjRice_et^ 
l 2013h . I 6 l Allers fc Liu 


Although G1 866 ABC is a triple system, iRoias-Avala et al.l (|2012r) find it to be 
close to Solar metallicity ([M/H]=0.05±0.12) with a temperature similar to single 
M5 dwarfs according to analysis of an unresolved R~2700 iF-band spectrum. 
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Table 2. Spectral Fitting Results’^ 


Spectral 


T.// (K) 



log(g) 


Type 

YJH 

YJ 

H 

YJH 

YJ 

H 

Field Objects 

Ml 

3820tt“^ 

38661^13 

3738131^ 

^■^^-0.83 

A c;c;+0.98 
+■1+1-1.00 

4-581197 

M5 

2969;|i 

29341?^ 

34041“^ 

4 76+0-87 
+■ ‘O-1.06 

4 ‘15+1-14 
+"+l_o.99 

4-4ilS:^^ 

M 8 

266i;“^ 

263911^9 

29191313 

4 70+0-85 

+■ 'O'_0.93 

4 41++07 
+■+-‘--0.98 

4 80+1-19 
+-'1U_o,94 

LI 

23051^9 

231711“ 

2433119 

r no+0-53 
^‘*-*^-0.98 

4-551};°^ 

4 73+0.91 
+• ‘'1-1.23 

L5 

209814^ 

215911 

222911“ 

r . or+0.22 

^-'J^-0.38 

4.98l?;i 

4.94115? 

L 8 

20261^^ 

2187136^ 

21491^^? 

c 00+O.I6 
o.oo_0 21 

5.3811^3 

c 1 1+0.40 

T2 

16631342 

1374lJ?^ 

17821^31 

c 9-] +0.20 
^-^-‘--0.21 

c 4 c:+0.33 
'1-+'1-0.41 

5-20l9:?3 

T4.5 

121511 

12641119 

13281J38 

r 04+0.21 
^-'4+-0.23 

4 01+0.45 
+■11-‘--0.52 

5-4819;“ 

Young Objects 

M 8.57 

22461^7 

2337 I 9 I]® 

2347 I 57 ? 

5 70+0-18 
I*-'11-0.35 

4 ‘lO'’'1-^i 

+■'+1-0.94 

4 34+1-18 
+-'l+-1.00 

LO 7 

20931^^ 

22051^^° 

22391^^^ 

r Qq+0.07 

4.03li;33 

4.741?;?? 

L 57 

178011 

19691J1? 

1634137 

r 20+0-21 
o.^y_0 27 

5-0ll?ii 

4.801?;?? 

L7 VL-G 

1739131 

19471“ 

165911 

5 02+0-^0 
ll-ll^-0.30 

5.4419;“ 

4.65l?;?7 


^■6681 fit parameters are 50% quantile values, and quoted uncertainties are 
derived from 16% and 84% quantiles, equivalent to l-ci uncertainties for Gaussian 
posterior distributions (see Section 1 ^ . 
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Fig. 1.— Observed (dashed lines) and simul ated (red bars) neardnfrared spectra of fi eld M, L, and 


T dwarfs from IRTF/SpeX library spectra (jCushing et al. 


2005 


Ravner et al.ll2009l ). The binned 


spectra simulate the wavelength coverage and spectral sampling of P1640 YJH spectra (Phase 2). 
The calculation of the simulated spectra is described in Section YI?2\ The spectral template sources 
are listed in Table [Ij 
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Fig. 2.— Observed (dashed lines ) and simulated ( red bars) near-infrared spectra of young M 
and L spectral type objects (left; mice et al.l l2nind : iKirkpatrick et al.l 120061: iFahertv et al.l 120131 : 


Allers fc Lhjl2013l : iLiu et al.ll2013l ) and their closest field-age counterparts in our sample (right). 


The binned spectra simulate the wavelength coverage and spectral sampling of P1640 V JH spectra. 
The calculation of the simulated spectra is described in Section [221 The spectral template sources 
are listed in Tabled! 
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Fig. 3.— Synthetic spectra calculated with PHOENIX model atmosphere code trimmed and binned to 
match P1640 spectra as described in Section [3l The top panel shows spectra for the complete range 
of model temperatures in 100 K increments at a fixed surface gravity log(g)=5.0 [cgs] expected for 
field objects. The dusty models are shown for a fixed ^e// >2000 K, and cond models are shown 
for Teff <1900 K, hence the gap between spectra apparent at ~1.4 /rm. The middle panel shows 
spectra for the range of surface gravities at 2100 K (approximately late-M/early-L spectral types) 
from dusty models, and the bottom panel shows the same at a fixed re//=1200 K (approximately 
early-mid T dwarfs) from cond models. 
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Fig. 4.— Goodness of fit values for model spectra as a function of effective temperature (left) 
and surface gravity (right) for the simulated P1640 YJH spectra from the M5 (top), LI (middle), 
and T4.5 (bottom) spectral templates. This represents the first step in the fitting process and 
provides an overview of the complete parameter space. In these examples the lowest goodness of 
fit values are for model parameters 2950 K and logg=5.3 [cgs] for the M5 template, 2300 K and 
logg=5.5 [cgs] for the LI and 1200 K and logg=5.3 [cgs] for the T4.5. These values are used as 
starting points for the MCMC procedure.The overview of parameter space provided by this step in 
the htting procedure demonstrates that the low resolution near-infrared spectra of very low mass 
objects are sensitive to temperature (as expected) and are increasingly sensitive to surface gravity 
for later spectral types. 
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Fig. 5.— MCMC results for the field M5 spectral template object for: complete YJH spectrum 
(purple), YJ spectrum (blue) and H spectrum (red). The lower-left panel shows 1-, 2-, and 3-a 
contours of the posterior distribution functions for both parameters, while the top and right panels 
shows the histograms for temperature and surface gravity, respectively, marginalized over the other 
parameter. The distributions for the YJH and Y J spectra are much narrower in temperature than 
in surface gravity. The YJH and Y J-band spectra produce similar temperature results, but the 
i^-band distribution has an only slightly sloping tail toward higher temperatures. 
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Fig. 6.— Same as Figure [5] for the field LI spectral template object. The results for the complete 
VJH spectrum are better constrained in temperature than in surface gravity. Temperature is 
symmetric, while gravity has a longer tail to lower values. Y J-band spectra result in a more 
precise temperature estimate, while Lf-band and the complete YJH spectrum provide a more 
reliable gravity estimate as indicated by the sharp decrease in probability at gravities higher than 
~5.6, which are unphysically high but allowed in the model grid, and the less pronounced decrease 
in probability for log(g)<5.0. 
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Fig. 7.— Same as Figure [5] for the field L8 spectral template object. The distribution for the 
Y JH spectrum is narrow and symmetric in both temperature and gravity. The Y Jdistribution 
has a broader, but still single-peaked, distribution in temperature but increases toward unphysical 
high values in surface temperature. The 77-band distribution is broad in temperature, covering 
the entire range of both the YJH and YJ distributions in temperature. In gravity the i7-band 
distribution is consistent with the YJH distribution with a long, flat tail toward lower values. 
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Fig. 8.— Same as Figure [5] for the field T4.5 spectral template object. The results from different 
spectral bands are generally consistent, with long tails to higher effective temperatures for the 
individual V J- and //-band spectra and to lower surface gravities for the Y J-band spectrum. 
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Fig. 9.— Best fit PHOENIX model spectra for simulated P1640 spectra (gray error bars) of field 
M, L, and T dwarfs. Vertical error bars represent the constant noise value with SNR=10 at the 
peak flux, and horizontal errors bars represent the width of one wavelength channel. Purple lines 
represent the best fit model spectrum using all flux points in the simulated spectra except four 
points in the H 2 O band 1.34-1.42 |Um). Blue dashed lines represent the best fit model spectra 
using just the ~yj-band flux points (0.995-1.31 /rm) and red dot-dashed lines, the 77-band flux 
points (1.45-1.80 /rm). M and L dwarf fits use the dusty version of the PHOENIX models, and T 
dwarf fits use the cond version. 
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Fig. 10.— Best fit PHOENIX dusty model spectra for simulated P1640 spectra (gray error bars) of 
young M and L spectral type objects (left) and their closest field-age counterparts in our sample 
(right). Vertical error bars represent the constant noise value with SNR=10 at the peak flux, and 
horizontal errors bars represent the width of one wavelength channel. Purple lines represent the 
best fit model spectrum using all flux points in the simulated spectra except four points in the H 2 O 
band 1.34-1.42 //m). Blue dashed lines represent the best fit model spectra using just the ~VJ- 
band flux points (0.995-1.31 urn) and red dot-dashed lines, the i/-band flux points (1.45-1.80 
M and L dwarf fits use the dusty version of the PHOENIX models. 
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Fig. 11.— Best fit results in effective temperature versus spectral type for simulated P1640 data 
fit as yj-band (blue upward triangles), i7-band (red downward triangles) and YJH-haiid (purple 
squares) spectra. Field objects are represented by filled symbols and young objects, by open 
symbols. The gray filled reg i on sh ows temperatures derived for obje c ts wi th the same spectral 
types from iLuhman fc: Riekd (jl999l l for M dwarfs and IStenhens et al.l (|2009l l for L and T dwarfs 
with a range of ±100 K. 
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Fig. 12.— Best fit results in surface gravity versus spectral type for simulated PI640 data fit as Y J- 
band (blue upward triangles), i^-band (red downward triangles) and YJH-hand (purple squares) 
spectra. Field objects are represented by filled symbols and young objects, by open symbols. The 
solid shaded region sh ows the range of su rface gravities for f ield objects (age > 5 00 M yr) predicted by 


evolut ionary models: ISiess et al.l (j2000l ) for T^f f >3000 K.lBaraffe et al.l (120021 ) and IChabrier et al 
(|200nl i [DUSTYOO] for 3000>Te//>1500 K, and iBaraffe et abl (l200,3l l jcOND03] for Te//<1500 K. 
The hatched gray regions shows predictions from the DUSTYOO models for effective temperatures 
~2500-1500 K, which corresponds approximately to the spectral types of young objects in our 
sample. The surface gravity results for field objects are closer to predictions from evolutionary 
models and more consistent between individual wavelength regions than the results for the young 
objects. 
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Fig. 13.— MCMC results for young and field objects of approximately the same spectral type, ~M8 
objects on the left and L5 objects on the right, field objects on top and young objects on the bottom. 
Colors correspond to different spectral fits as in Figure O The peak of the distributions in effective 
temperature for the young objects are cooler than the corresponding peak for the field objects. 
Both young objects also have distributions that peak at higher surface gravities than predicted by 
evolutionary models using ages constrained by their likely membership in nearby young moving 
groups. See Section ISTl for a complete discussion of the results. 
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Fig. 14.— Same as Figure[5]for the young L7 VL-G object. The distributions from different spectral 
regions are substantially less consistent than for the field objects. For example, the 2-cr contours 
of the Y J-band (blue) and i7-band (red) distributions just barely overlap. While the i7-band 
distribution peaks at a surface gravity predicted for a young, late-type object (see Figure [12]), the 
Y JH and Y J-band distributions are peaked at higher surface gravities and have only a long tail 
toward lower values. 
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Fig. 15.— Best fit effective temperatures (top) and surface gravities (bottom) for fits to the 
simulated P1640 spectra of the M5 (squares), LI (diamonds), and T4.5 (triangles) dwarf templates 
with varied SNRs, as described in Section 15.71 Symbols on the bottom panel are offset slightly 
in log(SNR) for clarity. The error bars show the standard deviation of the distribution of best fit 
parameters for 10^ resampled spectra using comparisons with the non-interpolated model grid. The 
results suggest that SNR>5 is required for temperature determination and SNR~10 is optimal for 
inferring surface gravity from best fit model spectra. 























































